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ABSTRACT 

Aims. We aim to demonstrate that the Herschel ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying 
candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals 
a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our 
approach by comparing the properties of the most likely candidate disks to those of the known population. 

Methods. We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio 
method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance 
cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS A'-band images. 
Results. We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely 
debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties 
to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search 
sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and 
outline the future prospects for our debris disk search programme. 
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1. Introduction 
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! The Herschel ATLAS or H-ATLAS (lEales et all 1201 Oh is 

• the largest Open Time Key Project on the Herschel Space 
Observatory (Pilbratt et al., 2010), and will ultimately map over 

. 500 square degr e es with both the PACS an d SPIRE instruments 
; (lPoghtschetal.L 120101: iGriffinetaLl l20Ta) . H-ATLAS is de- 
I signed to revolutionise our view of dust and dust-obscured star 

formation by detecting ~ 250000 galaxies in the far-infrared. 

The primary goal of the H-ATLAS is to study galaxy forma- 
, tion and evolution (see the other H-ATLAS articles in this vol- 

• ume), however the unrivalled sensitivity and wide-area coverage 
' means that H-ATLAS can also reveal dust in a range of more 

local (i.e. within the Milky Way) environments. At the sensitiv- 
ity hmits of H-ATLAS (a Scr threshold of 33 i riJy at 250 ^um 
', measured from the scie nce demonstration data, IPasquale et al.l 
1 1201 tf , 'Rigbv et al]|2010 ) this implies the potential to detect ana- 
logues of the welLknown debris disks (e.g. Holland et al., 1998; 
iGreaves et all 1 19981) out to distances of between 20 and 150 pc. 
A search for debris disks in H-ATLAS offers a powerful 
complement to those deeper and more targeted studies that are 
cuiTently being undertaken with Herschel (DUNES, DEBRIS 
and GASPS — see publications in th is volume), are set to 
be carried out with SCUBA-2 (SUNS: iMatthews et"^ l2007h 
and have been perfomed by Spitzer (thoroughly described in 
[Carpenter et all 12009, and references therein). With its wide 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA 



areal coverage H-ATLAS is a shallower tier to these studies, 
but a tier with the potential to search a much larger number of 
stars for bright debris disks and with a large body of s upport- 
ing high quality optical and infrared legacy data (Eale s et aU 
2010). T he H-ATLAS fi elds are covered by SDSS DR7 in ugriz 
( Abazaii an et al.L l2009l) and U KIDSS Large Area Survey in 
YJHK ( Lawrence et aI.L l2007l) with forthcoming deeper INT 
and VST KIDS optical data, VISTA VIKING in the near- 
infrared, WISE in the mid-infrared and GMRT radio contin- 
uum. The supporting optical & infrared data allows straightfor- 
ward selection of main sequence star s in the H -ATLAS fields 
via the main sequence colour locus (Covev et al., 2007) and the 
use of techniques to exclude contaminating background galax- 
ies, such as the inspection of deep /T-band imag es for extended 
objects and the use of the FIR-radio correlation dCarilli & YunL 
1999). Finally, a debris disk search in H-ATLAS is completely 
serendipitous and carried out in parallel with the primary science 
programme. 

In Sect|2]we show that the full H-ATLAS survey will allow 
us to search ~ 10000 main sequence stars for the presence of 
bright debris disks analogous to Beta Pictoris and ~ 1 000 stars 
for Fomalhaut analogues. This large search sample means that 
H-ATLAS will be much more sensitive to rarities in the debris 
disk population than targeted surveys , leading to st ringent tests 
of stochastic disk evolution models (IWvattL l2008h and poten- 
tially uncovering bright and/or co ld disks that may have under- 
gone recent disruptive events (e.g. lLisse et al.Ll2007h . H-ATLAS 
will allow us to answer questions such as how frequent are bright 
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systems such as Beta Pictoris or HR 4796 and is there an upper 
limit to the amount of debris formed during disk evolution? 



2. Photometric selection of main sequence stars in 
the l-l-ATLAS Science Demonstration Field 

The H-ATLAS science demonstration (SD) field occupies 
roughly 16 square degrees and is centred at 09:05:30.0 
+00:30:00 (J2000). Descriptions of the PACS and SPIRE im- 
ages obtained in par allel mode and th e da ta reduction procedure 
used can be found in'Ib ar et all (1201 Ol) and lPasquale et al.l (1201 Ol) 
respectively. Point sources were identified within the images us- 
ing a combin ation of PSF filtering, Gaussian fitting and aperture 
photometry (Rig bv et al.L l2010l) . The median 5cr noise values 
(including confusion noise) in the SPIRE 250, 350 & 500 ;um 
maps are 33, 36 & 45 mJy/beam respectively. The PACS images 
are roughly a factor 2 noiser than predicted with Scr noise levels 
of 105 and 138 mJy/beam at 100 andl60^ m respectively. The 
H-ATLAS source catalogue for the SD field ("Rigbv et all 1201 Ol) 
contains 6878 band-merged sources with flux density > 5cr in 
at least one of the 5 H-ATLAS bands (100, 160, 250, 350, 500 
fim). The catalogue excludes higher noise regions at the edge of 
the map and thus covers an effective area of 14.5 square degrees. 
Each source in the H-ATLAS catalogue has been matched to the 
corresponding most likely SDSS DR7 object within a 10" search 
radiu s. This matching process is described further in lSmith et alj 
(1201 0) and is an implementation of the Bayesi an Likelihood 
Ratio technique of lSutherland & SaundersI (Il992h . 

In order to search for debris disks in this catalogue, we 
must first identify a sample of main seque nce stars. We use th e 
main sequence colour locus identified by ICovev et aP (l2007h . 
which constrains the location of main sequence stars in SDSS 
colour space. This approach allows us to maximise our search 
sample by going to faint magnitudes and takes advantage of 
the well-calibrated and well-understood SDSS optical photome- 
try. We_use_a_4;dirnensional main sequence locus as described 
in iKimball et alJ (l2009l) rather than the full 7-dimensional 
SDSS+2MASS locus of lCovevet al] (|2007|) . The infrared ex- 
cess of a warm debris disk at /T-band can move our target stars 
away from the nominal locus, and as our aim is to identify stars 
that are potential debris disk hosts, we thus do not use 2MASS 
or UKIDSS colours in our photometric selection. 

Figure [T] shows a colour-magnitude diagram of the 180000 
star like objects (selected with 'probPSF^l') detected by SDSS 
DR7 in the H-ATLAS SD field. The general population are 
shown by red dots and those falling within 2 'u nits' of the 4- 
dimensional main sequence color locus defined bvKimb all et alJ 
(|2009|) are shown in blue. Note that the colour locus is 4- 
dimensional and Fig. [T] shows only a ID cut through the locus. 
Figure [U shows that the bulk of the main sequence stellar popu- 
lation detected by SDSS is comprised of faint and likely distant 
halo stars. However, there are a substantial number of relatively 
bright (11 < i < 17) and hence likely nearby stars for which it is 
possible that H-ATLAS could detect associated debris disks. 

We estimate the maximum distance to which H-ATLAS 
could be sensitive to debris disks by scaling from the Spitzer 
MIPS and SCUBA SEDs of known examples (Be ta Pictoris, 
iRebu lletal .1120081: Epsilon Eridani. iBackman et a lJ 12009; HR 
8799, Lisse et al. 2007; F omalh aut,[Stapelfeldt et alJ|l2004l Vega 
& HR 4796, Sheretetani200 4'). The most sensitive wavelength 
of H-ATLAS is 250 fim (at which the stellar photospheric con- 
tribution is minimal) and we find that the maximum distances 
for these debris disk analogues are: HR 4796 <200 pc. Beta 
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Fig. 1. SDSS / vs. / - z colour magnitude diagram of SDSS point 
sources (i.e. with probPSF = 1) in the H-ATLAS SD field. A 
magnitude cut of / < 21.3 has been applied to exclude sources 
with large photometric error. Dots indicate field SDSS point 
sources (red) and point sources within the main sequence stel- 
lar locus (blue). The grey shaded polygon indicates the colour- 
magnitude region occupied by main sequence s tars between dis- 
tances of 4 and 200 pc (calculated using the iDavenport et al.l 
(2006) photometric parallax relation). Stars within the main se- 
quence locus that fall into this box are shown as larger green 
points for clarity. 



Pictoris <150 pc, Fomalhaut <80 pc, Vega/HR 8799 <50 pc and 
Epsilon Eridani <20 pc. We indicate the colour-magnitude re- 
gion in which main sequence dwarf stars at a distance of 4- 
200 pc should lie in Fig. [T] by a grey shaded box, calculated 
using the distance modulus and the M, vs. (; - z) p hotometric 
parallax relation for dwarf stars from IDavenport et al. (2006!). 
Note that late M and L dwarfs follow a shallower relation for 
(/ - z) > 1 .26 (IWest et al.l, 120051), and also that we have extrap- 
olated the Davenport et al. (2006) relation to (/ - z) < 0.2 to 
account for the bluer stars in the sample. 

In the H-ATLAS SD field for maximum distances of 200, 
150, 80 and 50 pc (i.e. sensitive to HR 4796, Beta Pictoris, 
Fomalhaut and Vega analogues) we find a total of 851, 340, 31 
and 9 stars respectively on the main sequence locus. There are 
no stars on this locus nearer than 20 pc in the SD field, which is 
likely an effect of the SDSS becoming saturated for near, bright 
stars. Such stars can be obtained from the Tycho-2 & Hipparcos 
catalogues (lH0g et all l2000t Ivan LeeuwenL l2007h . although in 
this paper we focus upon the larger and better selected SDSS 
sample of main sequence stars. Assuming similar stellar densi- 
ties in the remaining 550 square degrees that will be mapped in 
the full survey, H-ATLAS will thus encompass a search sam- 
ple on the order of 10000 main sequence stars for the brightest 
debris disks (HR 4796 and Beta Pictoris analogues) and on the 
order of 300-1000 stars for Fomalhaut and Vega analogues. 



3. Candidate debris disks in the H-ATLAS SD Field 

We identify candidate debris disks by taking a sample of H- 
ATLAS sources from the catalogue that have a high reliabil- 
ity match (>80% reliability) to SDSS DR7 catalogued sources 
(Smith et al., 2010). There are 2334 sources that pass this crite- 
rion. We then filter this list to only include SDSS point sources 
CprobPSF^l') and identify point sources on the 4-dimensional 
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Fig. 2. SDSS / vs. / - z colour magnitude diagram of SDSS point 
sources with high reliability matches to H-ATLAS catalogue 
sources. Objects whose SDSS colours place them on the stellar 
locus are identified by open circles. Thre grey shaded region in- 
dicates the colour-magnitude region occupied by stars between 4 
and 200 pc as in Fig.[T] Solid and dashed lines indicate distances 
of 300 & 400 pc respectively. The two candidate debris disks are 
identified by their H-ATLAS catalogue ID. 



main sequence colour locus described in the previous section. 
We find a total of 204 H-ATLAS sources matched to SDSS 
sources with 'probPSF=r, of which 78 sources fall within the 
main sequence locus (see Fig. |2]i. Within this sample we ex- 
pect considerable contamination from dust obscured QSOs or 
unresolved galaxies whose optical colours are reddened into 
the main sequence locus (Ivezic et al., 2002i). Indeed, the sam- 
ple have a median r-band magnitude of 1 9.8, fainter than th e 
value at which galaxies dominate over stars (ICovev et al.L 120071) . 
Unsurprisingly due to their optical colours, none of the H- 
ATLAS main sequence locus objects have measured SDSS spec- 
troscopic redshifts which would allow us to select out QSOs and 
unresolved galaxies. The inclusion of UKIDSS ne ar-IR colours 
in the selection would aid in this discrimination (llvezic et al.L 
I2OO2I) . but as we mention in the previous section, would also 
select against possible /T-band excess in debris disk stars. 

To identify candidate disks in our sample we apply a pho- 
tometric distance cut to select the brightest and nearest objects 
that are least likely to be QSOs or unresolved galaxies and the 
most likely to be debris disks. We apply an initial photometric 
distance cut of 200 pc to select the most likely candidate disks, 
though Fig. |2]shows that there are a further 7 objects at 200^00 
pc distance that could be massive or luminous disk candidates. 
The candidates that pass the main sequence colour locus and 
photometric distance tests are then finally subject to a detailed 
inspection of SDSS DR7 ugriz and UKIDSS LAS YJHK images 
to reject the presence of possible contaminating galaxies, and the 
wider field of the H-ATLAS SPIRE images to search for contam- 
inating cirrus. We stress that our search technique reveals candi- 
date disks. Spectroscopic confirmation of the host star spectral 
types, more accurate spectrophotometric distances, higher res- 
olution PACS or SCUBA-2/ALMA imaging, or scattered light 
imaging are required to confirm these objects as debris disks and 
to better constrain their physical properties such as temperature 
and mass. 

We focus our following discussion on the two closest candi- 
date disks found within a photometric distance cut of 200 pc: H- 



ATLAS J0903 15.8+015758 and H-ATLAS J090240.2-014351. 
As we will show, the physical properties of these objects are 
within the spectrum of known debris disks and the H-ATLAS 
detections are thus consistent with a debris disk hypothesis. 
We summarise the FIR and stellar properties of H-ATLAS 
J090315.8 and H-ATLAS J090240.2 in Tableland present 
three colour images of the disks in Fig. [3] The g - i colours 
of the host stars imp ly spectral types of K2 and G5 respectively 
dCovev et all l2007h . Their 2MASS colours or brightnesses are 
inconsistent with those of giant stars (J - H & H - Ks for both 
stars is <0.3 and < 0.1 respectively). The photospheric flux of 
these stars at 250 fixn is of the order of a few fj.}y and so we 
can be confident that the 250 /urn emission is a genuine excess 
over the stellar spectrum. Both disks are unresolved at 250 //m, 
although H-ATLAS J090315.8 shows a marginal extension to 
the West. H-ATLAS J090240.2 on the other hand is compact 
and centred on the star's position to within the pointing accuracy 
of Herschel. Inspection of UKIDSS /T-band images shows that 
background galaxy contamination is unlikely. 

Both candidate disks are detected in the 250 & 350 //m 
SPIRE bands, but not at 100, 160 or 500 //m (both disks have 
emission significant at the 2cr level at 500 yum). With only two 
flux points it is difficult to constrain either the fractional lumi- 
nosity or the temperature of the disk. We estimate temperature 
by fitting a fixed yS=0.5 modified blackbody (e.g. I Wvatt et all 
|2005) to the 110, 160 //m upper limits and 250, 350 & 500 
//m flux measurements. As the 110 and 160 //m flux are only 
an upper limit the derived temperatures should be considered 
as strict upper limits. We note that as our p hotometric dis- 
tance estimates are only good to within ~50% (Davenport et alj, 
2006), the error in derived disk mass is largely dominated by 
distance ra ther than temperature . Us ing the standard tec hniques 
outlined in [Holland eTal] ( Il998l) and lSheret etal] (|2004 we de- 
rive masses for the two candidate disks of 0.06-0.13 M^ and 
0.5-1. 3Me for H-ATLAS J090315.8 and H-ATLAS J090240.2 
respectively. The mass of both candidate disks is within t he ob- 
served range of k nown disks (IWvattl 120081 : iLisse et aU 120071 : 
ISheret et aU, l2004t) . comparing specifically to HD 12167 with a 
mass of 1 Me and j6 Pic with a mass of 0.1 Me. dSheret et all 
12004 

4. Conclusions & future prospects for H-ATLAS 
debris disk searches 

We have described a search method for debris disks in the 
Herschel ATLAS and present the two most likely candidate disks 
in the H-ATLAS SD field: H-ATLAS J090315.8 and H-ATLAS 
J090240.2, ~ 0.1 and 1 Me mass candidate disks around KO and 
G5 stars respectively. We also identify a further population of 
76 SDSS point sources that are associated with FIR/sub-mm 
emission and whose optical colours place them on the main se- 
quence locus. The majority of these objects are likely to be dust - 
obscured QSOs or unresolved galaxies (e.g. llvezic et al.Lr2002h . 
although the brighter and nearer (7 objects lie within a photo- 
metric distance of 400 pc) may be potentially luminous debris 
disks. Follow-up spectroscopy is required to investigate these 
hypotheses. What is clear from the H-ATLAS SD field is that 
bright disks are rare — we have searched 85 1 stars within 200 
pc for FIR/sub-mm emission and find only two candidate disks 
brighter than 33 mJy at 250 fim. Our search finds a much lower 
fraction of candidate debris disks than previous Spitz er and 
SCUBA studies (e.g. ICarpen teret al.. 2009; Hillenbran d et all 
2008; Wyatt, 2008), which typically find a 7-14% disk detec- 
tion fraction. However we note that direct comparisons between 
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Table 1. F IR and stellar prop erties of the candidate disks & their associated sta rs. Spectral types are est imated from the (g - i) colour 
relation of lCovev et alJ (12007 ) & distances from the photometric parallax from'Paven port et aP (12006 '). Temperatures are estimated 
by a gr eybody fit to the 25 0, 350 & 500 yum fluxes and the upper limits of 140 mJy at 160 /im, and 105 mJy at 110 jim with a fixed 
(5- 0.5 (jWyatt et al., 2005). Masses are calculated using the standard technique with a mass absorption coefficient k - 1.7 cm^ g"^ 
at 850 yum (.Sheretet al.>,2004,) . 



H- ATLAS ID 



SDSS ID 



Sp. type Distance (pc) Fiso (miy) F350 (miy) F500 (miy) T (K) Dust mass {Mq) 



J0903 15.8+015758 
1090240.2-014351 



J090315.91+015800.9 


K2 


90-130 


51 ±14 


53 ± 14 


19 ±9 


65 


0.06-0.13 


1090240.10-014349.7 


G5 


190-290 


86 + 12 


48 ±7 


17 ±9 


60 


0.5-1.3 




Fig. 3. Three colour images of the two candidate debris disks 
detected by H-ATLAS. Red & green channels are 2MASS J & 
Ks respectively. The blue channel is SPIRE ISOj-im, smoothed 
by a 3 pixel Gaussian kernel to increase the signal to noise ratio. 
Contours are of 250 fim emission starting at 3cr and spaced by 
2cr. 



these detection fractions should not be made as H-ATLAS is 
flux-limited rather than volume-limited and we do not reach the 
exq uisite photospheric sign al-to-noise levels of the Spitzer stud- 
ies dCarpenter et all l2009h . Hence H-ATLAS is likely to only 
detect the bright end of the debris disk population. 

The search that we present in this Letter is a forerunner to 
a much more extensive programme that will be carried out in 
the future. The full H-ATLAS will contain a search sample of 
~ 10 000 photometrically selected main sequence stars out to 200 
pc, allowing us to place more stringent limits upon the frequency 
of the bright end of the debris disk population. With a large and 
well-selected sample of main sequence stars covering a range 
of spectral types we will also be able to carry out a stacking 



analysis (e.g. lKurczvnski & GawiserLl2010l) to determine the sta- 
tistical frequency of disk occurrence as a function of spectral 
type. We also plan to include bright stars from the Tycho-2 & 
Hipparcos catalogues, which will extend our search to nearer 
main sequence stars. A preliminary search of the Tycho-2 cat- 
alogue indicates that none of the 569 stars found within the H- 
ATLAS SD field are associated with detectable debris disks. For 
the 7 nearest Hipparcos main sequence stars within our field 
(which lie between 30 & 80 pc) this means that H-ATLAS can 
place upper limits on their disk masses of 0.01-0.07 M® (or be- 
tween 0.8 and 5.7 Lunar masses of dust, assuming Tdust = 40 K). 
Finally, H-ATLAS has signifi cant legacy potential for the GAIA 
mission (Lindegren et al., 2008), which will determine distances 
and spectral types for all stars brighter than r = 20 in the H- 
ATLAS fields. In combination with GAIA parallaxes, H-ATLAS 
will be able to determine precise disk mass upper limits for a 
large sample of stars. 
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